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I describe a new likelihood technique, based on counts-in-cells statis¬ 
tics, that I use to analyze repeating in the BATSE IB and 2B cata¬ 
logues. Using the IB data, I find that repeating is preferred over non¬ 
repeating by 4.3:1 odds, with a well-defined peak at 5-6 repetitions per 
source. I find that the post-lB data are consistent with the repeating 
model inferred from the IB data, after taking into account the lower 
fraction of bursts with well-determined positions. Combining the two 
data sets, I find that the odds favoring repeating over non-repeating 
are almost unaffected at 4:1, with a narrower peak at 5 repetitions per 
source. I conclude that the data sets are consistent both with each 
other and with repeating, and that for these data sets the odds favor 
repeating. 
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INTRODUCTION 


By studying clustering in the angular distribution of bursts, Quashnock and 
jmb (1,2) found evidence that 7 -ray burst sources repeat. Using a nearest 
ftfeghbour analysis of the BATSE IB catalogue (3), we found a significant 
teJtcess of neighbours on angular scales less than 5°. Since this is less than 
median positional error of bursts in the catalogue, we concluded that 
jWrst sources repeat multiple times, on a time-scale of months. Wang and 
■^ijigenfelter (4) presented evidence that five particular bursts arise from a 
rSgle repeating source: They also found evidence of repeating by studying 
sJSmial and temporal correlations in the IB data (5). Using a particular model 
(©repeating, Strohmayer, Fenimore and Miralles ( 6 ) also found evidence for 
Ideating 

^^learly it was paramount to extend these analyses to the BATSE 2B cat¬ 
alogue (7). Unfortunately, the failure of the tape recorders on board the 
Compton Observatory led to a decrease of ~ 1/3 in the fraction of bursts with 
well-determined (so-called “non-MAXBC”) positions in the new (or 2B-1B) 
data. While two-point angular correlation function and nearest neighbour 
analyses of the 2B catalogue failed to confirm repeating ( 8 ), this is expected 
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given the positional errors of order 7° and the aforementioned drop in ex¬ 
posure (5). This is particularly the case if the typical number of observed 
repetitions from each source is small (9). Assuming that the bursts in the 
IB catalogue form a fair sample, Lamb and Quashnock (10) simulated what 
would be expected in the 2B-1B and 2B catalogues. We found that repeat¬ 
ing is not detectable, due to the lower fraction of bursts with non-MAXBC 
positions. 

If burst sources do repeat, this would significantly constrain the range of 
allowed models and favor a Galactic origin. Given these implications and the 
controversy surrounding the issue of repeating, it is clearly imperative to use 
the most powerful and sensitive techniques to test the repeating hypothesis - 
in particular because of the drop of exposure in the 2B catalogue. Likelihood 
methods (11) are known to be the most sensitive and give the best possible 
determination of model parameters. 

Here I describe a new likelihood method based on counts-in-cells that I 
have developed and used to analyze the clustering of gamma-ray bursts in 
the BATSE IB and 2B catalogues (9). It allows us to test various repeating 
models parametrized by the number Ar of repeating sources and the number 
Arep of bursts emitted by each repeating source, and to include the important 
effects of exposure and positional errors. 

COUNTS-IN-CELLS LIKELIHOOD 

Let Aceii be a large number of cells, each of fixed solid angle size fi, each 
centered on a random position on the sky. Let Cjv to be the number of 
these cells having A bursts in them, where A = 0, 1, 2,... . I then define the 
observed counts-in-cells distribution, 

Pn = Cj^fNcell, ( 1 ) 

which is the probability that a randomly chosen cell of size Q has A bursts 
in it. 

The counts-in-cells distribution contains information about clustering on 
scales comparable to the angular size of the cell. Indeed, the expected dis¬ 
tribution Qjv is Poisson when the bursts are uniformly distributed on the 
sky: 


= (2) 

where Af is the number density of bursts. 

I have calculated the expected counts-in-cell distribution Qjv for a repeating 
model specified by the parameters and Arep, where is the number 
of repeating sources each of which bursts exactly Arep times. By definition 
Arep > 2. Note that this is the actwaf number of repetitions per source, not the 
observed number, which is on average much smaller because the BATSE sky 
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exposure e is significantly smaller than unity [e = 0.34 for the IB catalogue 
and only 0.26 for the 2B-1B (non-MAXBC) catalogue]. I include both the 
effect of finite positional accuracy (0err = 6.8° is the median error) and of 
sky exposure when calculating the expected counts-in-cells distribution in the 
repeating model. I also allow for an integer number N^t of background sources 
that only burst once. 

Once I calculate the expected counts-in-cells distribution Qjv for a given 
set of parameters of the model, I ask how likely is the observed distribution 
Pjv for the data, given such a distribution. This likelihood C is given by the 
following (9): 


log £ = Aceii ^ Pn log Qn + const.. 

N 


(3) 


I then use Bayes’ Theorem to interpret the likelihood in terms of a probability 
distribution for the integer values of the model parameters; namely the number 
of repeating sources and the number of repetitions for each repeating 
source. 

I have analyzed the BATSE IB and 2B catalogues using this counts-in-cells 
technique and a cell size of 5°. Figure 1 (left panel) shows the 1-, 2- and 3-cr 
contours in the (Arep, Ar)-plane. The cross marks the maximum likelihood 
location. While the credible interval for the number of repeating sources is 
broad, that for the number of repetitions is considerably narrower. This is 
shown in Figure 1 (right panel), where I have marginalized over the number 
of repeating sources. There is a well-defined peak at 5 or 6 repetitions per 
source. The repeating model is favored over the non-repeating model by odds 
of 4.3:1. 

Figure 2 (left panel) shows the results of the same analysis for the 2B-1B 
(non-MAXBC) catalogue. Note that the contours in the (Arep, Ar)-plane are 
now much larger and extend down to smaller numbers of repeating sources 
and repetitions per source. Indeed, the odds of repeating versus non-repeating 
have fallen to 0.85:1 (basically equal odds). Nevertheless, the 1-cr credible re¬ 
gions in Figures 1 and 2 largely overlap, and the maximum likelihood locations 
are quite close. Figure 2 (right panel) again shows the probability distribution 
for the number of repetitions per source. While the peak has now fallen to 2 
repetitions, the probability of 5 repetitions is almost as large. 

I have also combined the two data sets, and find the results shown in Figure 
3. Interestingly, the contours of the 1-, 2-, and 3-cr credible regions for the 
2B catalogue (left panel) are actually smaller than for the IB. This is evident 
in the right panel, which shows a well-defined peak in probability at 5 repeti¬ 
tions per source. Combining the two data sets, I find that the odds favoring 
repeating over non-repeating are almost unaffected, at 4:1, relative to that 
found for the IB. 
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FIG. 1. (left panel) 1-, 2- and 3 -it contours in the (JMrep, i\^r)-plane from analysis 
of the IB catalogue, (right panel) Probability distribution of iVrep from the same 
analysis. 




FIG. 2. S ame as Figure 1, from analysis of the 2B-1B (non-MAXBC) catalogue. 
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FIG. 3. S ame as Figure 1, from analysis of the combined catalogues. 


CONCLUSIONS 

Considering all of these results, I conclude that the first and second data 
sets are consistent both with each other and with repeating, and that for 
the combined data sets the odds still favor repeating. Thus the 2B-1B data 
alone neither confirm nor refute the repeating hypothesis, nor can they given 
the drop in exposure relative to the IB catalogue (5). The 3B-2B catalogue 
containing 570 new bursts is not expected to suffer from this, and it should 
offer a fair test of the repeating hypothesis. 
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